Order from disorder in the double exchange model 
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The magnetic excitations of the double exchange (DE) model are usually discussed in terms of 
an equivalent ferromagnetic Heisenberg model. However this equivalence is valid only at a quasi- 
classical level — we show that both quantum and thermal corrections to the magnetic properties 
of DE model differ from any effective Heisenberg model because its spin excitations interact only 
indirectly, through the exchange of charge fluctuations. We also find that the competition between 
ferromagnetic double exchange and an antiferromagnetic superexchange provides a new example of 
an "order from disorder" phenomenon — an intermediate spin configuration (either a canted or a 
spiral state) is selected by quantum and/or thermal fluctuations. 



•4-* 



C 

O 

> 
O 

OS 

m 



o 
o 

-i— > 



T3 

c 

o 
o 



X 



Many magnetic systems of current experimental inter- 
est, for example the colossal magnetoresistance (CMR) 
Manganites and Pyrochlores || , consist of itinerant 
electrons interacting with an array of localized magnetic 
moments with spin S. The simplest models of these sys- 
tems comprise a single tight binding band of electrons 
interacting with localized core spins by a ferromagnetic 
(Hund's rule) exchange interaction Jh t 



Hi 



E 

(ij)a 



C ia C ja 



Jh ^ S, 

ia/3 



(1) 



where the sum (ij) is restricted to neighboring sites. At 
a classical level, the ground state of this model in D > 1 
must be ferromagnetic, since the Hund's rule requires 
the spin of itinerant electrons to be locally aligned with 
the core spins, and the kinetic energy of electrons is in 
turn minimized by making all the electron spins parallel. 
This effect is usually called "double exchange" , and in 
this context the model of Eq. ([!]) is referred to as the 
double exchange ferromagnet (DEFM). 

At a quasi-classical level the spin dynamics in a DEFM 
can be described by an effective nearest-neighbor Heisen- 
berg model —J\ Y](jj\ Si-Sj with ferromagnetic exchange 
integral J\ = t/AS 2 , where t is the expectation value of 
the kinetic energy per bond in the lattice [Q-|],[l2| . fn ad- 
dition, core spins also interact via a direct superexchange 
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In CMR materials, J2 is believed to be positive (antifer- 
romagnetic), and so competes with the DE mechanism. 

A natural question to ask is how a system described by 
TL = Tii + 7^2 evolves towards an antiferromagnet with 
increasing J 2 ? As observed by de Gennes M , for classical 
spins (S — > 00), the FM becomes unstable at J1/J2 = 1. 
At larger J2, the competition between J2 and the ki- 
netic energy gives rise to an intermediate phase where 



the neighboring lattice spins are misaligned by an angle 
9, where cos 9/2 = J1/J2 (in this regime, the DE model 
is no longer equivalent to Heisenberg model) . However, 9 
alone does not specify a particular intermediate configu- 
ration — there exist an infinite set of classically degener- 
ate states with the same 9. The two ends of this set are 
the two-sublattice canted phase and the spiral phase [Q . 

We performed a spin-wave analysis for the canted 
phase of the DEFM and found ui sw (q) = for all q. This 
result reflects the local symmetry linking all the differ- 
ent classical groundstates — we can take any (classical) 
spin on the A sublattice of the canted phase and rotate it 
about the direction of magnetization of the B sublattice 
without changing the angle 9. Since it costs no energy to 
make such an excitation, the system cannot distinguish 
between different states, and is magnetically disordered 
even at T = Q [Q. This argument, however, does not 
hold at finite S, and we anticipate that quantum and/or 
thermal fluctuations will enable the system to choose its 
true groundstate. Such "order from disorder" effects have 
been widely discussed in the context of magnetic insula- 
tors, e. g. for Kagome antiferromagnets |l3| — but not 
in the context of the metallic DEFM. 

The most direct route to an answer to what kind of 
order is preferred in the DE model is to determine the 
first instability of the DEFM, by looking for the wave 
vector Q* at which the spin-wave dispersion becomes 
unstable. The quasi-classical arguments leading to an 
effective Heisenberg model cannot answer this question 
since they predict a FM spinwave spectrum of the form 
Wsw(q) = 2zS(Ji — J2XI — 7 g ), where z is the lattice coor- 
dination number, and 7 9 = (1/z) e tkS where {6} runs 
over nearest neighbors. Thus uj sw (q) vanishes identically 
for J2 — Ji, and no special wave vector is singled out. In 
order to find the true groundstate for J2 > J\, we then 
must go beyond the usual mapping of the DEFM onto 
an effective Heisenberg model. 
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In this paper we argue that the equivalence between 
DEFM and the nearest-neighbor Heisenberg ferromag- 
net holds only for 5 — > oo, while fluctuation corrections 
in the DEFM are governed by fermions and are different 
from those in the Heisenberg model. We show analyti- 
cally that for both small densities of electrons (x -C 1) 
and small densities of holes (1 - i « 1), the first insta- 
bility of the DEFM is against a two-sublattice canted 
structure with Q* — (tt,tt, 7t). At intermediate densities 
(0.3 < x < 0.9 in 3D), the first instability is against a 
spiral spin configuration with Q* w (0,0,0). Similar re- 
sults hold for the 2D case. We also show that thermal 
corrections to uj sw (q) compete with the quantum correc- 
tions, and the trade-off between these two effects gives 
rise to re-entrant transition between canted and spiral 
phases with increasing T. 

The noncquivalence between DEFM and the Heisen- 
berg ferromagnet has been earlier detected in numerical 
studies ■ Analytically, this has been demonstrated 

at T = by Golosov |§. Our T = results agree with 
his, but we present more physical understanding. Our 
calculation scheme and results at T ^ are entirely new. 

We now turn to the calculations. Given that the the 
Hund's rule coupling Jjj is the largest energy scale in 
the problem it is desirable to diagonalize this term first 
and project out all electrons not locally aligned with the 
core spins. We should also define spin-wave excitations 
such that they are the true Goldstone modes of the or- 
der parameter, i. e. transverse fluctuations of the com- 
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posite spin S l = S t + \ Y, a 

plish both of these goals by the following procedure : (i) 
introducing new Fermi operators fa and pi which cre- 
ate electronic "up" and "down" states aligned with the 
quantization axis of the composite spin, and (ii) general- 
izing the Holstein-Primakoff transformation to the case 
where the length of the spin is itself an operator, intro- 
ducing a corresponding bosonic operator a. We work in 
the large 5 limit usual for spin-wave theories. Full details 
of the transformation will be given elsewhere here 
we simply state that the transformation and its inverse 
are unitary, and satisfy all required (anti-)commutation 
relations e. g. {/,/ f } = {p,p f } = [a,a f ] = 1, [/>»] = 
{fiP} = e tc We have also checked that our expansion 
procedure reproduces all features of an exact solution of 
the DE model on two sites Jig] , 

Although the full analysis of the order from disorder 
phenomenon required us to derive the transformation to 
order 1/S 2 , to capture the essential physics it is sufficient 
to know the inverse transformation to order 1/5. To this 
order we have 

a = a I ! + V}v_p_ \_ Pp^ + 0(1/s3/2) 
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FIG. 1. a)-d) Diagrams contributing to spin-wave self en- 
ergy in Ferromagnet to order 0(1/S 2 ). Only diagrams a) and 
b) are physically relevant. e)-f) The representation of the di- 
agram b) via an effective four-boson interaction mediated by 
the charge suscpetibility (particle-hole polarization bubble). 
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where a is the Holstein-Primakoff boson associated with 
the core spin. Substituting this transformation into (Q), 
we find that the Hund rule term in ([[]) does not contain 
any spin operators and reduces to 
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Clearly, p operators describe high-energy excitations and 
can be safely dropped. Simultaneously, the hopping term 
in ([!]) transforms into the Hamiltionian which describes a 
single band of spinless fermions interacting with (initially 
dispersionless) Holstein-Primakoff bosons: 
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The perturbation theory for the bosonic propagator 
D(q,fl) is straightforward. We have D~ 1 (q, f2) = 17 — 
S(g, f2). The self-energy physically comes from the fact 
that the oscillations of the core spins destroy perfect 
alignment of electron spins, and this increases electron 
kinetic energy. To 0(1/5), £((?, f2) comes from a single 
loop of fermions Fig. la and evaluates to 
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FIG. 2. Schematic corrections to the spin-wave dispersion 
in the (1, 1, 1) direction at zero and finite temperatures and 
small electron density x. Dotted lines show zero temperature 
classical (0(1/ S)) dispersion; dashed lines show renormalized 
classical dispersion at 0(1/ S 2 ); solid lines show spectrum in- 
cluding all quantum/thermal effects at 0(1/S 2 ). Insets show 
residual dispersion at instability against canted (T = 0) and 
spiral (T 7^ 0) phases 



so at this level we indeed reproduce the spin-wave spec- 
trum of a nearest-neighbor Heisenberg FM. 

The calculation of the spin-wave dispersion at higher 
order requires some care. There are four self energy con- 
tributions at 0(1/S 2 ), one from each of the diagrams in 
Fig. la-d. At finite temperature these include a contribu- 
tion from a six-fold term in 7i t , omitted in Eqn. |^. How- 
ever there is considerable cancellation between terms, 
and all important physical effects come from the dia- 
gram Fig. lb, which can be thought of as an effective 
four-boson interaction mediated by Fermions (Fig. le- 
f). Assembling all contributions to the self-energy, split- 
ting it into quantum and thermal pieces, and neglecting 
the frequency dependence of the self-energy (which is 
0(1/ S 3 ) ), we obtain a renormalized spin-wave disper- 
sion fi(g) = 2zJiS(l - 7,) + E^L (g) + S^ 2) ( 9 ) where 
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tively. The latter can be approximated by step functions 
as we will be interested in T < J\S -C t. 

We now analyze the form of T,^ 2 \q), beginning with 
the case T = 0. The first term in (|J) simply renor- 
malizes the classical Heisenberg-like spinwave disper- 
sion, while the second has a dependence on q which is 
quite different from that in the nearest-neighbor Heisen- 
berg model. This term is either positive or negative 
throughout the Brillouin zone, depending on the elec- 
tronic density, and is symmetric under q — ► 7? — q where 
7? = (tt, 7T, 7t). Along the zone diagonal, it reduces to 
&I(x) [1 - cos(2 9 )] where I(x) = (l/N 2 ) n F (p)(l- 

«f(0) hp - 7p +7r / 2 )/(7 P - li) changes sign twice as a 
function of x. This form of the correction to Heisen- 
berg dispersion is comparable to that found in numerical 
studies of the DE model on a ring S . 

We find that for intermediate densities 0.31(7) < x < 
0.92(7) in 2D, and 0.31(7) < x < 0.94(2) in 3D, 
I(x) < and quantum effects cause a relative soften- 
ing of spinwave modes near the zone center. This means 
that the first instability of the DEFM with competing 
AF exchange interactions is against a spiral phase with 
Q* w (0, 0, 0). On the other hand, for a small density of 
electrons (or holes), I(x) > 0, and quantum effects in- 
stead lead to a relative softening of modes near the zone 
boundary. In this case the spin-wave spectrum first be- 
comes unstable against a canted spin configuration with 
Q* = (7r,7r,7r). These results are in perfect agreement 
with earlier studies jll| . Finally, we note that if we for- 
mally extend our large S analysis to arbitrarily small S, 
there exists a critical value of spin S — S* « 1 (in 3D) for 
which the DEFM becomes unstable even in the absence 
of AF interactions. This opens up the possibility that 
for small S (e. g. , S = 1/2), the ground state of the DE 
model may not be a FM, as suggested by some numerical 
studies JTot - 

We now proceed to finite T. In a Heisenberg model, 
finite temperature effects do not change the form of spin- 
wave dispersion, but the overall scale of the dispersion is 
reduced by a factor cx T 5 / 2 /S 2 [jL6|. In the case of the 
DEFM, the result is more complex and coincides with 
the behavior of a Heisenberg FM only in a very lim- 
ited range of frequencies and temperatures. Indeed, there 
are three typical momenta in our problem — an exter- 
nal q, a fermionic p F , and a typical bosonic momentum 
hyp ~ {T / JiS) 1 / 2 (in units where the lattice constant 
a = 1). Consider for defmiteness the case of small x 
when pf = ( l OTr 2 x) 1 /' i <C 1. Evaluating the self-energy 
(^) by expanding all 7^ factors to second power in mo- 
menta, we obtain in 3D, for q -C p F , (T/J1S) 1 / 2 



x(7p-m/2 - 7p+?-m/a)(7p+m/2 ~ Ip+q-m/i)- 
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Here rh = q — I, and ns(q) — ns(0(g)) and np(k) = 
nF(efe) are Bose and Fermi distribution functions, respec- 



where $(^) is a smooth function with <^(v — > 0) — 
1, $(v = 1) = 0, $(1/ > 1) « -1/V 2 , and 
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FIG. 3. Phase diagram showing first instability as a func- 
tion of carrier concentration x and temperature T. For some 
range of dopings x\ < x < X2 a re-entrant phase transition is 
predicted. 
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(2) 

see that Sf ; (g) agrees with the Heisenberg model only 
at (T '/ ' JiS) 1 / 2 < pf- In the opposite limit, wc have 



£f J (?) oc q 2 pp(T/J 1 Sf/ 2 . Moreover, the sign of T,^>(q) 
changes between the two limits. 

For large external q: q^> pp, (T/JiS) 1 ^ 2 , we obtain 
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We see that at these q, the temperature correction 
to the spin-wave energy is positive, and is of order 
pp{T/ J\S) b l 2 . It depends only weakly on q at interme- 
diate momenta, and vanishes at q — (ir, tt, tt). This result 
is indeed very different from that in the Heisenberg FM. 
The temperature correction to the spin-wave dispersion 
is shown schematically in Fig.^Ja. A very similar behav- 
ior holds in 2D case, the only difference is that the 2D 
analog of the scaling function $(f) does not change sign 
at v ~ 0(1) but rather behaves as l/v A at large v. 

The different nature of quantum/thermal effects in the 
DEFM and the Heisenberg FM is a result of the differ- 
ent nature of the coupling between spins. In a Heisen- 
berg FM, the q 2 T b / 2 form of the temperature correc- 
tion comes from the fact that four-boson vertex scales as 
q 2 l 2 In the DEFM, the interaction between spin 

waves is mediated by the charge susceptibility of the 
Fermi gas, U(q - I) (see Fig. (le)). At q, l typ <C pf, 
n(<7 — I) can be approximated by n(0), and the effective 
interaction has the same form as in the Heisenberg FM. 
However, when \q — l\ 3> pf, the susceptibility decreases 
as n(m) tx p 2 F /m 2 . This replaces either one power of 
T ' I J\S or q 2 factor by p 2 F in the renormalization of the 
dispersion, exactly as we found. 

(2) 

The unusual temperature dependence of T, y T (q) gives 
rise to the possibility of a re-entrant transition between 
spiral and canted states with varying T. As discussed 
above, at small p F , quantum fluctuations favor canted 
phase. On the other hand, classical fluctuations at low 



{T/JiS) 1 / 2 < pf soften the dispersion near q = and 
hence favor the spiral state. As a result, the transition 
line between the two states bends towards smaller x at 
finite T. However, as T increases and becomes larger 
than JiSp F , the sign of (q) changes, and thermal 
fluctuations now favor the canted phase. As a result, the 
transition line now bends towards higher x with increas- 
ing T. This can give rise to a re-entrant transition — 
when T increases at at a given, small x, the canted state 
first becomes unstable towards the spiral state, and then 
returns back at even larger T. A possible phase diagram 
for x <C 1 is shown schematically in Fig.||b. The region 
occupied by the canted phase first shrinks and then ex- 
pands with increasing T. Qualitatively similar behavior 
will hold for 1 - x < 1. 

To summarize, in this paper we introduce a novel large 
S expansion scheme for systems with strong Hund's rule 
coupeling, and use it to show how and why the spin dy- 
namics of a DEFM differ from those of a Heisenberg FM. 
We find that the two models are equivalent at the classi- 
cal level (S — > oo), but that both quantum and thermal 
corrections in the DEFM are different because the inter- 
action between spinwaves is mediated by fermions. We 
also find that the DEFM provides a new example of an 
"order from disorder" phenomenon — in the case of a 
competition between ferromagnetic double exchange and 
a direct antiferromagnetic supcrexchange, the classical 
ground state is infinitely degenerate, but fluctuations lift 
the degeneracy and select a true intermediate spin config- 
uration, either a canted state or a spiral state, depending 
on the electron density. This together with the calculated 
temperature corrections gives rise to the possibility of an 
unusual re-entrant transition between spiral and canted 
states with varying T . 

The softening of zone boundary spin-waves which we 
predict has been observed in neutron scattering experi- 



ments on the CMR Manganites |I7|, but at a hole doping 
x = 1 — x ~ 0.3 for which our zero temperature theory 
would predict a relative hardening of the dispersion at 
(■7T, 7r,7r). It has been suggested that this softening may 
be due to the influence of optical phonons , or orbital 
degrees of freedom [|l9|. We note, however, that the DE 
mechanism does predict a softening at the physically rele- 
vant doping, for high enough temperatures. A careful ex- 
perimental examination of the temperature dependence 
of the spin-wave spectrum is therefore necessary to de- 
termine which mechanism contains the relevant physics. 

The issue left for further studies is a possible phase sep- 
aration in the non- ferromagnetic regime [pp| . To study 
this possibility in our approach, one has to analyze the 
sign of the longitudinal susceptibility in, e.g., canted 
phase. If it is negative, then the system is unstable to- 
wards phase separation pl[ . There calculations are cur- 
rently under way. 
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